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1. DEVELOPMENT OP A MATHEMATICAL MODEL FOR STEADY STATE ISOELECTRIC 
FOCUSING. 

Isoelectric focusing (IEP) is a high resolution electrophoretic 
technique which is widely employed for protein separations. The technique 
is based upon the observation that proteins exhibit zero mobility at 
unique positions (their isoelectric points) In a pH gradient. The intro- 
duction by LKB (under the tradename Ampholine) of mixtures of syntheti- 
cally random polyaminopolycarboxylic acids provided an easy means of 
establishing a broad pH gradient and contributed to a large extent to 
the remarkable popularity of the analytical technique. Preparative 
IEE has enjoyed less popularity as a result of (a) scaling constraints 
of the classical focusing apparatus, (b) the loss during fractionation 
of the resolution achieved during focusing, and (c) the unavoidable 
contamination by Ampholine of the fractionated proteins. We are 
addressing the last of these problems and proposing the elimination of 
Ampholine from the system by establishing the pH gradient with simple 
ampholytes. To that end we have developed a mathematical model of 
simple ampholyte IEF. A simple ampholyte, in our terminology is one 
whose behavior near its isoelectric point results from the influence 
of only two ionizeable groups. Predictions of the model and computer 
simulation of focusing experiments are of direct importance to our 
efforts to obtain Ampholine-free, stable pH gradients. 

1.1. Simple Model based on the Electroneutrality Approximation. 

Our first effort to develop a model of IEF was largely 
directed by Svensson’s pioneering work‘d in the early 1960 f s. He 
considered the focusing of a single ampholyte component, and we general- 
ized the treatment to two, three and multi- component systems. This 

^H. Svensson, Acta Chem. Scand. 15 , 325-341 (1961) 
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approach employs the common form of Ohm* s law to describe the electric 
field as a function of the solution conductance. The model consists 
of an ordinary differential equation for each ampholyte which is 
a statement of conservation of mass, and a coupling algebraic equation 
which represents the electroneutralxty approximation. The model 
equations were solved using discretization techniques based on the 
Runge-Kutta-Merson integration rule and the solution was implemented 
on a digital computer using DARE P simulation software , The results 
of the solution were presented as the pH, conductance, and the concen- 
tration of each ampholyte as a function of distance along the column 
axis. This model was used for the computer simulations discussed 
later in this report. The salient aspects of the model are described 
in a paper, "Theory and Computer Simulation of Isoelectric Focusing 
with Simple Ampholytes", submitted to the LKE Symposium — Electrofocus 
78 and included in Appendix A. The developemnt of the multi-component 
system is included as Appendix B. 

As a result of comments from Dr. Dudley Saville during a visit 

3 

to our laboratory and subsequent communications , we have reappraised 
this model. Dr. Saville indicated that the contribution of the current 
due to diffusion should be considered in calculating the effective 
electric field, i.e., a more complete expression of Ohm's law must he 
employed for a rigorous treatment. The effect of this correction 
will of course be dependent upon the magnitude of the diffusion current, 
relative to the applied current. Our preliminary calculations indicate 
that this is not large for the cases considered; however, we feel it 
is necessary to specifically include the diffusion current in the 

2 

J. J. Lucas and J, V. Wait, DARE P - A Portable CSSL-Type 
Simulation Language, Simulation (1975) 

3 

D. A. Saville, personal communication 
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general form of the model and neglect its contribution only when it 
can be shown to be negligible. 

1.2. Revised General Model of Steady State Isoelectric Focusing. 

The revised model describing a one-dimensional IEF system 
is based upon the following assumptions; 

(1) the system is at steady state, i.e. there is no net 
flux of any ampholyte component; 

(2) the concentrations of component subspecies are described 
by equations of chemical equilibria, i.e. chemical reactions (proton 
associations and dissociations) are rapid relative to transport processes; 

(3) the system is isothermal and there is no macroscopic flow. 
1.2.1. General Equations. 

Flux (F^) of i-th species 

-F. - e w.z.n.V0 + RT w.Vn. (1) 

x x x x r x x 

w. - mobility (m^sec - ^V 

1 3 
n^ - concentration (moles/m ) 

z^ — valence (dimensionless) 

0 - local electric potential (V) 

R - universal gas constant (8.314 joules /°K . mole) 

4 

e - molar charge (6.9366.10 Q) 

T - temperature (°K) 

Q - Coulomb's unit change 

Conservation of mass (conservation of i-th ions) 

V-F.jiSL - 0 (2) 

- rate of production of i-th species per unit 

volume (rate of change of charge density) O^r -) 

3 sec 
m 

Poisson’s equation 

2 N 

e O « -e 2 z.n. C3) 

.-.xx 

X=1 

e - dielectric constant (5.48 •lO - ’’^ ) 

V.nr 

N - number of species 


3 “ 


Conservation of charge (net production of charge is zero) 


N 

2 z.R. - 0 
xi 


C4) 


Conservation of j-th ampholyte 


1= 3 2 

2 R. ™ . 0 

i~jV 1 


(5) 


R R 

j, 5 j . +1 . . . . , R . - rate of production 

xx j 2 

of species j^+1, . . , , , j^) - component subspecies 

of j-th ampholyte 
1.2.2. Model Construction. 

The conservation of charge equation: Multiplying equation 

(2) by the valence and adding these equations for all species gives 


N 


N 


2 z.V.F . + 2 z.R. 

• t X X ...XI 
1=1 1=1 


( 6 ) 


Taking into account (4) the above equation becomes 
N 


2 z.V.F. = 0 
x x 

i=l 


(7) 


or 


N 

V. 2 z.F. 

.-XX 

i=l 


= 0 


(3) 


Since valences z. are constant. 

i 

From here we will treat only the one dimensional case, 

when the symbol V represents the derivative, that means we substitute 

V = — 
dx’ 

Equation (8) can be integrated once, giving 
U 

2 z.F. - -J 


i=l 


1 x 


J - total current density 
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Using (1) in the above we get 

N 9 ji* . N dn. 

e(2 i.Vi.jffi + SI S z.w. T • 

1=1 1 11 & 1=1 " xdx 


Ampholyte conservation: A procedure similar to one used 

in above can be applied here. The equations (2) for component sub- 


species of j-th ampholyte are added, which results in 


i=j 2 dF i=j 

2 y — ~ + .2 R. = 0 
dx . x 


(10) 


Using (5) in the above equation leads to the following relation 
d i== 0 2 

dx , S . *1 - 0 OH 

1=3 -j_ 


The equation (10) can be integrated, which results in 


a.-3 9 ■ 

2 F. = c, 
1 ; 


3 = '1,2, .... ,M 
M - number of ampholytes 
The flux of any ampholyte is zero 


c. = 0 
3 


since we assume steady state. 


This results in the following M equations 


2 F. - 0 
. . x 

1=3 1 


Using (1) in (14) gives 


1 3 2 dd i“3 7 dn, 

Z z w n ) % + ET 2 2 w. -=-3 
*_4 xxx ax . . x dx 
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The equation (9) and equations (14) constitute a system of JHi 
equations with W unknowns (n^yn^, .,,n^). For model completeness, 
some other relations are needed since N" > IfKL. 

Chemical Equilibria: 

We assume that the chemical equilibria are obtained rapidly with 
respect to transport processes, and that the concentrations of com- 
ponent subspecies are related by the dissociation constants. To 
illustrate this approach let us consider an ampholyte A^ which 
dissociates into simole valence ions i.e. At and A. . 

■ 33 


We use the following notation 
.+ 


Species : H ’ OH 

Concentrations: n- n, 


• • * # 


1 2 *** 


At 

3 


n» n. n. 
3 1 3 1 "®"1 3 < 


A. + H 
3 


A. + H 
3 



K. • = 
J 1 

4 

— “ A. ; 

c 2 3 

C 3 

C 3«fc rt A- 

II 

♦ 

^2 


The dissociation constants are 
n-n. 

K. 1 


n. 


+1 


K. 


J 2 ^ 


n-n. 

1 3- 


n. 


C 16 ) 


K = n.n„ 
w 12 


It is possible to use (16) in order to express n. _,n. and n n in 

J l +1 J 2 

terms of n- and n. , thus reducing the number of variables in the 
3 1 

model. We can proceed similarly in the cases of more complex ampholytes, 
which dissociating lead to more then 3 component subspecies. This 
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way we can. express the concentrations of component subspecies of 
an ampholyte in terms of the hydrogen ion concentration (n^) and 
either the concentration of one chosen component subspecies or 
the total ampholyte concentration (n^ ) 

j 

n 

t. = I n. 
i . . i 


( 17 ) 


1.2.3,' Isoelectric Focusing Model With a Single Simple Ampholyte 
Example. 

Notation 


I: 

Species : 
Concentration: 


H 


2 

QH“ 

n„ 


3 

A 

n„ 


A 


n, 


lw. , 

H_0 — H + OH - 
2 c 2 

. + fi. 

A + H *>* — A 

C 4 


K = n„n 
w 


12 


K 2 = ^ 


A + H + A + 

c 6 


n„n 


K. = 


13 

n, 


= 


n l n 5 

n„ 


K l = 


(18) 


K 


w 


Z 1 = 1 » z 2 


-1, z 3 = 0, ^ 


n l n 2 

1, z 5 = -1 


For this example we have only one ampholyte (M=l) and 3 

From the equation (9) we obtain 

^2 dfl ^ dn_ 

■ eC J z i *±«i> s + “ l vi $r ■ 


-J 


(19) 


* jo= 5. 


( 20 ) 


Similarly from (15) we have 


5 j fi 5 dn- 

s(S 5 +HI» l3r 


= 0 


( 21 ) 


Using (IS) and (19) we can write 


K 


n„ = 


w 


2 n„ 


Ur “ 


ti^K^ 


(22) 


n, = 


n l n 3 


4 ^ 

From (22) we compute 


dn r 

i 

dx 


K dn- 
w 1 

2 dx 
n- 


dn, n„ dn, 
4 _ _3 1 

dx dx 


n^ dn. 




(23) 


dn c K„ dn_ K« dn,, 
Z n^ i _j z J 

dx 2 dx n- dx 

JL 


Using (22) and (23) in (20) we obtain after some simple 


manipulations 

K 


n. 


K 2, in 3 


_3 ■ x K 2°3 ' ^ ,74 

tWi+Wj 2 + » 4 + w- 5 2 J dx + ( K_ V w s n 5 dx 

n^ l 1 1 


_ -J e , ' K w , n l n 3 n 3 K 2 v d0 

RT “ RT ^ W l n l + w 2 n x + W 4 K ± + w 5 ti ± ^ dx 


( 24 ) 


Similarly using (22) and (23) in (21) we obtain 

w^ n^ dn^ w^ d- 

( k 7 n 3 " W 5 K 2 2 5 dx“ + tw 3 + kT *1 + W 5 nT^ ± 

In, 1.1 


■ 8 - 


(25) 


e 


RT 


, . n l n 3 „ n 3 x 
< w 4 “ * 5 ^ 5 ^ 3 


d0 

dx 


Equation 3 becomes 


, e, v d 2 0 

(- /e) 1 - ^ - n 2 + n z 

dx 


n c 


(26) 


Taking into account the value^ of the constant ( £ /e = 6.06 x 10 ^) 

and assuming the rate of change of the electric field is much smaller 
15 -2 

than 10 V m , we replace Poisson 7 s equation by an algebraic one, 
representing the electroneutrality approximation. Using (22) we write 

n_ - K " 4n.ii, - K 9 (n-/nL.) « 0 (27) 

1 w n^ 1 3 2 3 1 


1.3. Direction of Continuing Research. 

We have a mathematical model of steady state IEF based 
on the electroneutrality approximation which has been written in 
the forms for single and multicomponent systems. We have a revised 
general mathematical model of steady state XEE, which has been written 
for the single ampholyte system. We intend to solve the equations 
for the single component system in both forms and compare their 
solutions and results. Using the most advantageous approach, we 
will then address multiple component focusing systems. Wo major 
difficulties are forseen at the level of two and three components. At 
this level we will attempt quantitative experimental verification of 
the model and conduct remodeling if required. The feasibility of 
achieving a solution with several ampholyte components will be evalu- 
ted and implemented if considered reasonable. We intend to express the 


4 

Cells 


S. R. Vaccaro and H. 
(manuscript) 


S. Green, Ionic Processes in Excitable 
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equations in terras of the total concentration of each ampholyte. 

We also intend to evaluate the feasibility of extending the 
model to consider the transient state. A transient state solution 
would, in addition to enlarging the theoretical foundation of iso- 
electric focusing, provide additional insight into the mechanics of 
the approach to steady state. This may allow further exploitation of 
the process as a kinetic rather than static event, and lead to increases 
in time efficiency. This last factor may be quite important as the 
time required for focusing the gradient can be considerable. Dr, 

D. Saville will be providing us with collaborative assistance during 
the coming year, 

2. COMPUTER SIMULATIONS OF STEADY STATE IEF 

2.1* Simulations and Their Evaluation. 

Computer simulations are intended to serve three functions 
in our research. First, the simulations increase our understanding 
of the focusing process by allowing us to vary ampholyte characteristics 
or experimental conditions and observe their affect on the IEF system. 
Second, the simulations provide theoretical data against which the 
experimental data may be compared to establish the validity of the 
theoretical model and its parameters. Third, the method of simulation 
will hopefully be developed to the level that it may provide the 

basis for designing focusing systems with specifically desired 

! 

characteristics . 

Computer simulations have been conducted with a variety of ampholyte 
systems at the level of two and three components using the simple model. 
The simulations are set up by entering specific values for the applied 
current, the initial (boundary ) concentrations of the ampholytes, and the 


physicochemical characteristics for each ampholyte which define its elec- 
trophoretic behavior, i.e., its proton dissociation constants, mobility, 
and diffusion constant. In some cases these parameters have been 
given values which reflect a real compound, while in other cases values 
representing purely imaginary components have been used. The latter 
method allows sequential variation of a given paramenter (or set) to 
establish what effect this variation will have on the characteristics 
of the focusing system. The index of the library of computer simulations 
we have conducted is included as Appendix C. For the comparison of 
related systems, the boundary concentrations of each component may be 
established as different (to represent a limit of the column) or as 
equal (to represent some intermediate position) . When using equal 
concentrations of both components as the boundary condition in a two 
component system, the simulation may be directed toward either the 
anode or the cathode. To represent a given experimental setup, the 
cross-sectional area of the column and its length must be specified 
to allow integration of the ampholyte concentrations over these dimen- 
sions. The integrated values are then the total amount of ampholyte 
contained within the column volume. The simulations have routin el y 
been conducted using the dimensions of the ISCO Column Model 212 which 
is one of the focusing apparatus routinely employed in our laboratory. 

The results of a simulation using two components representing 
glutamic acid and histidine are illustrated by the concentration 
profiles (figure 1) and the conductivity and pH profiles (figure 2) . 

The conductivity and pH profiles are of course of primary interest in 
regard to utilization of the focusing system, while the concentration 
profiles are of interest in that they represent the manner in which the 
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CONCENTRATION CmM> 



•t 


COLL 


Fig.. 1: The concentration p 
using glutamic acid (solid I 
current enroloyed was 4 x 10“ 



CONDUCTANCE Crnhos/crn 


* 

s- 


0 

t 


COLUMN AXIS (cm) 

Fig. 2: The conductivity (solid' line) and pH 

(broken line) profiles of the two-component IEF 
system. The conditions were identical to Fig. 1 


pH and conductivity are. established. The profiles obtained from a 
series of simulations may be replotted together to demonstrate the 
effect of parameter variation. The pH and conductivity data presented 
in figures 3 and 4 were generated by setting the pK T s of one ampholyte 
(pKg^= 6 and pK^- 8) while varying the pK T s of the other ampholyte 
systematically so that the A pK, i.e, pK^- pK^, constant 

while the pi . of the component varies. The boundary width, e.g. 
roughly the distance between points representing 5 % and 95% of one 
of the components maximum concentration, reflected in these profiles 
decreases at large values of pl^ - pl^ and increases for small values 
of this parameter. The pH data in figure 5 were generated by holding 
the pK f s of one ampholyte constant and varying the pK^ - pK^ value 
of the other without changing its pi. The boundary represented by 
these profiles sharpens at small values of pK^ - pK^ for either ampho- 
lyte and broadens for large values of this parameter . These figures 
also indicate that rather linear, shallow pH gradients may be established 
with ampholytes having either a small A pK (good ampholytes) or a 
large A pit (poor ampholytes) . The pH and conductivity data of figures 
6 and 7 respectively were generated varying the equal boundary con- 
centrations of both components. In this case the parameters were 
specified to represent one poor ampholyte (ApK=6, pi = 5) and 
one good ampholyte ( A pK - 2, pi = 7) . It can be seen that the 
higher concentrations produced more linear pH profiles over a narrow 
range and higher conductivity with a concomitant reduction in boundary 
sharpness. 

It was noted that the width of the boundary between two ampholyte 
components is inversely proportional to the amperage. The effect of 
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pi SERIES OPK CONSTANT) . pH 


•6/15 13s 43 pi “4 
• 6/15 13s 45 pi-5 


6/16 9* 47 P I -6 

6/16 9s S0 pi “6. 5 — 
6/16 9*51 pl**3 


Fig. 3: The pH profiles generated using one component with 

a pi of 7 , while the pK's of the other component were varied 
so that the pi ranged from 3 to 6.5. The current employed 
was 4 x 10"' A. 










ApK SERIES C P I CONSTANT!) pH 


-6/16 9:59ApK=8 
-6/16 l£L-00ApK=6 
•6/16 9: 57 ApK=4 
■6/16 g s S4ApK=2 


Fig. 5: The pH profiles generated by holding the pK’s of one 

component constant (pEC^-6 and pK_ 2 =8) , while the pK's of the 
other component (pl=5) were varied so that the value of pK^j 
ranged between 2 and 8. The current employed was 4x10"^ A. 






CM CO 




different amperages is demonstrated in figure 8, which gives the 
concentration profiles of two ampholytes representing glutamic acid 
and histidine from an equal concentration position. To extend the 
boundary over 20 cm requires a current of approximately 0.4 microamp. 

With a current of one milliamp the boundary width is reduced to 
approximately 0.1 mm. Thus at one milliamp the focused system is 
essentially an isoelectric zone of glutamic acid contiguous to an 
isoelectric zone of histidine through a very Sharp boundary expressed 
as a pH discontinuity. The sharp separation in this system more 
closely resembles an isatachophoretic separation than the type of sepa- 
ration commonly thought to occur in IEF . At comparable amperage when 
the intermediate component of a three-component system is present 
in excess of the two end components, an isoelectric zone in the middle 
of the column results. This isoelectric zone is expressed as bath a 
pH plateau and a conductivity gap. This phenomenon is illustrated 
graphically in figures 9 and 10. The concept of ampholyte concentration 
profiles resembling bell-shaped curves^ is clearly incorrect in this 
instance. A non-Gaussian ampholyte concentration profile and conduc- 
tivity gap have been observed experimentally in other laboratories 

5 

when modifying Ampholine gradients by amino acid addition . 

2.2. Direction of Continuing Research. 

Obviously since these simulations were conducted using the 
early model based on the electroneutrality approximation, the important 
aspects of the systems need to be reevaluated through simulations 
using the revised model. 

To date the simulated systems have been evaluated only qualitatively, 

5 

see for example, R. K. Brown, M. L. Caspers and S. N. Vinogradow, 
Carrier Ampholyte Distribution, in Electrofocusing and Isotachophoresis 
(eds. B. J. Radola and D, Graesslin) pp. 87-96, de Gruyter C1977} 
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CONCENTRATION CmM.) 


CURRENT SERIES, CONCENTRATION 


COLUMN AXIS (cm) 


Fig. 8: The concentration profiles of glutamic acid (CA) and 
histidine (CE) from an equal concentration boundary position. 
The current was varied from 0.4 to 3.6 microamps. 






CONCENTRATION CmM) 


COLUMN AXIS (cm) 


Fig. 9: The concentration profiles of the three-component 

system using glutamic acid {solid line) , histidine {dotted 
line), and a third component (broken line) with an intermediate 
pi (pK- =4, pK 0 = 6). The current employed was 1.7x10“6a. 


COLUMN AXIS (cm) 


Fig. 10: The conductivity (solid line) and pH (broken 

line) profiles of the three component system. The 
ampholytes and conditions are identical to Fig. 9. 





i.e. , with regard to conductivity, the best gradients are those 
exhibiting a relatively constant, high conductivity profile; with 
regard to pH, a linear profile would normally be preferred. It is 
expected that both narrow and broad range pH gradients are of interest; 
the narrow range (approximately one pH unit) would be preferred for 
separations within a specified region such as a preparative puri- 
fication of a single protein; the broad range 'pH gradient would be 
useful analytically or in preparative applications where high reso- 
lution is not required. 

To quantitatively evaluate focusing systems, one can resort to 
Rilbe's definition^ of the resolving power of the method: 


Apl 


D (d- 
= 3 E (- 

M 


(dpH/dx) 


- h 


-du/ dpH)j 

where Api is the pH difference between just-resolved protein zones, 
D is the diffusion coefficient of the protein and du/dpH is the 
slope of mobility vs pH at the component's isoelectric point. This 
could be applied by establishing a unit or constant value for 
(D/-du/dpH) and calculating maximum and mean values of Apl for 
each gradient. This of course assumes that the sample ampholytes 
are in low enough concentration that they have no influence on the 
established gradient. Alternatively, one of the ampholytes of the 
model stay be given the characteristics of a sample ampholyte, and 
its distribution in the focusing column evaluated. It would be 
advantageous to incorporate this type of analysis into the computer 
program.. 



Rilbe, 


Ann. N. Y. Acad. 


Sci. 209, 11-22 (1973) 
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3, EXPERIMENTAL VERIFICATION OF THE THEORETICAL MODEL 
3,1. Experimental Findings 

To confirm the theoretical model, we need accurate profiles 
of a system focused under several different experimental conditions. 

Since pH, conductivity, and ampholyte concentration profiles all 
represent different aspects of the same information, any of these may 
be utilized. The technique of analysis of course must be capable of 
discrimination within a small fraction of the length of the established 
gradient. We do not currently have an experimental setup which provides 
the requisite resolution, but we have made some qualitative observations* 
The most obvious testable prediction of the model is that the 
boundary between glutamic acid and histidine will be very narrow when 
the applied current is in the milliampere range. This was tested in 
two ways using the ISCO column: 1) the col umn was focused containing 

only the ampholytes glutamic acid (Glu, pi » 3.22) and histidine 
(His, pi = 7.58), then fractionated into 0.4 ml fractions (each 
representing 0.5 cm of column length) and the fractions analyzed; 2) 
the column was focused containing a visible marker in addition to the 
carrier ampholytes and isoelectric between them. 

In the former method, the fractions were analyzed as to conduc- 
tivity, amino acid analysis, absorption at 212 nm (histidine absorbance 
peak) , and pH. The absorbance and pH profiles of the column are pre- 
sented In figures 11 and 12. Although there is considerable loss of 
resolution during the fractionation procedure, the sharply discontinuous 
nature of the boundary is well-enough preserved to be illustrated by 
these data. 

When hemoglobin which is isoelectric around pH 7.2 was included 
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Fig. 11: The absorbance profile at 212 nm of the fractions eluted, 

from a focused column containing the ampholytes glutamic acid and 
histidine on a sorbitol density gradient. The applied current at 
steady state was 3 x 10”5 a. 
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Fig. 12: The pH profile of the fractions eluted from a focused column 

containing the ampholytes glutamic acid and histidine on a sorBitol 
density gradient. 


in the column as a visible marker, it focused in a very narrow band 
(less than 2mm wide) at the boundary between the Glu and his zones. 

When albumin dyed with bromophenol blue (which is isoelectric around 
pH 5,4) was included, both the hemoglobin and albumin bands were 
sharp and contiguous. When glycyl-glycine (pi = 5.65) was included 
as a third carrier ampholyte, the albumin focused between Glu and 
Gly-gly and the hemoglobin focused between Gly-gly and His in separate 
narrow bands. This demonstrates the possibility of separation of 
proteins by interposition or bracketing with carrier ampholytes 
analogous to mobility bracketing of proteins in isotachophoresis. 

When methyl red (pi 3.75) was used as a marker between His and Glu, 
it focused into a narrow band whose top half was red (the acidic 
color of methyl red) and whose bottom half was yellow (the dye’s 
basic color) . This strikingly demonstrated the very steep pH 
gradient occuring within the narrow boundary. While these, data 
are qualitative rather than quantitative, they illustrate the 
discontinuous nature of the boundary between adjacent ampholytes 
in this system. 

3.2. Direction of Continuing Research. 

To facilitate gathering the accurate profile data required 
to confirm the validity of the model, we are in the process of 
designing a system which will allow analysis along the length of 
a IEF column without requiring fractionation. Two types of apparatus 
are being considered: 1) a scanning conductivity probe which would 

move axially within the focusing column, and 2) a column capable 
of being scanned along its length in a spectrophotometer. The latter 
would utilize a small column contained within the lumen of a temperature 
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jacketed, quarts, flow-through, cell (Hellma 167-QS). We have 
conducted preliminary focusing experiments in the Hellma cell 
and achieved satisfactory results . 

4. IMPLEMENTATION OE SIMPLE AMPHOLYTE IEF FOR PROTEIN SEPARATIONS 

4.1. Results to Date. 

We have seen that simple ampholytes can be focused so that 
the boundaries between adjacent components are either broad, giving 
a pH gradient, or narrow, giving. a pH step. Broad boundaries are 
favored by low fields, high ampholyte concentrations, high values of 
pK^ -pK^ for biprotic ampholytes, and small differences in the iso- 
electric points of neighboring ampholytes, i.e. , pl^ -pl^. The 
converse of these conditions favors narrow boundaries. 

A pH-step system can be used to space or bracket specific 
proteins; however, the resolution of this technique is dependent 
upon the availability of ampholytes with closely spaced isoelectric 
point since an ampholyte of intermediate pi is required to separate 
sample proteins. To display a resolving power of 0.02 pH units 
would require 50 components with pi's spaced over a single pH unit. 

This sort of system is perhaps already available in the form of 
Ampholine, although the number of Ampbollne components used and 
their relative concentrations are of course not amenable to adjust- 
ment. Bracketing with high resolution using a reasonable number 
of components is probably more easily achievable in isotachophoresis 
than in IEF, since mobility is a flexible property for weak electrolytes 
and may be altered with pH, whereas isoelectric points are essentially 
inflexible. Bracketing with low resolution could conceivably be 
useful for separating a complex mixture into specific subgroups; 
however, it does not appear that such a technique would be of general 
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The use of pH-gradient systems is the alternative to using pH 
steps, and -we have begun to explore this alternative experimentally. 

In our preliminary attempts to focus proteins on simple ampholyte 
pH gradients we have primarily used a polyacrylamide slab gel 
system. The slab gels are easily manipulated and multiple gels can 
be run simultaneously. We have chosen to setup shallow, narrow- 
range pH gradients using two or three amino acids or dipeptides as 
carrier ampholytes. To date we have conducted experiments in only 
a few systems using the following amino acids and peptides: aspartic 

acid (pi = 2.77), glutamic acid (3.22), his tidy 1-glycine (6.81), 
histidyl-histidine (7.30), glycy 1-histidine (7.50), histidine (7.59), 
and 0-alany 1-histidine (8.17). The goal of these preliminary experi- 
ments was to determine if proteins can be focused with high resolution 
on a pH gradient formed from only two or three components. Por sample 
proteins we have used hemoglobin (Hb), catalase (Cat), or fractions 
of dialyzed egg protein (ED) . After focusing a strip of the gel is 
removed and sectioned, each section is eluted and the pH of the eluate 
recorded. The remainder of the gels are fixed with a Coomasie Blue 
solution to allow visualisation of the proteins . Protein distributions 
on representative gels and their pH profiles are presented as figures 
13 through 18. The pH range of the gel usually exceeds the pi range 
of the ampholytes, i.e. the. focusing column escapes ampholyte control 
at the ends. This phenomenon apparently is caused by the presence of 
buffers and salts In the system (by products of the acrylamide poly- 
merization reaction), and the magnitude of the effects is reduced 
by increasing ampholyte concentrations. While these results are very 
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GEL LENGTH v« pH 



PAG ISmM ASP* 16mM GLU* 100 Volte 
C P I»2„ 77) CpI-3. 22) 

WHOLE EGG WHITE PROTEIN 


12/4/78 

Fig. 13: The pH profile of an aspartic acid and glutamic acid 

gel with sample of egg white protein. 
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GEL LENGTH, cm 


GEL LENGTH vs pH 


ED7 



PAG 8mM HIS-HIS. SmM AI.A-HIS, 100 Volts 
<pl=7. 30) Cpl=8. 175 
HEMOGLOBIN 10mg/ml, EGG WHITE PROTEIN #7 

12/8/78 


Fig. 14: The pH profile of a his tidy 1-histidine and alanyl . 

histidine gel with samples of hemoglobin and egg white protein 
fraction. 



GEL LENGTH, cm 


GEL LENGTH vs pH 



pH 


PAG 8mM HIS-GLY, 8mM ALA-HXS, 200 Volta 
(pi -6. 81) C P I«8. 17) 

HEMOGLOBIN 10m S /ml, EGG WHITE PROTEIN #7 


Fig. 15: The pH profile of a histidyl-glycine and alanyl- 

histidine gel with samples of hemoglobin and egg white protein 
fraction. 
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GEL LENGTH ve pH 



P H 


PAG 8mM GLY-HIS, 8mM B ALA-HIS , 500 Volta 
(pi -6. 81) (pi “8. 17) 

HEMOGLOBIN lOng/ml. EGG WHITE PROTEIN 60mg/ral 
=2 hre. = 65 hre. 12/1/78 

Fig. 16: The pH profile of a glycyl-histidine and alanyl- 

histidine gel with samples of hemoglobin and egg white protein. 
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GEL LENGTH va pH 


Cat 



PAG 16mM HIS. lSmM ALA-HIS, 200 Volta 
(pi =7. 59) (pi =8. 175 

HEMOGLOBIN 10mg/ml. CATALASE 50mg/ml 

12/26/78 


Fig* 17: The pH profile of a histidine and alanyl -histidine gel 

with samples of hemoglobin and catalase. 



GEL LENGTH, cm 



PAG 16mM HIS-GLY, 16mM HIS, 16mM ALA-HIS 
<pl=*6. 81) C P I«7- 59) CpI-8. 17) 

HEMOGLOBIN 10m g /ml, CATALASE 50mg/ml, 200 Volte 


12/27/78 


Fig. 18: The pH profile of a his tidy 1-glycine, histidine, 

and alanyl-histidine gel with samples of hemoglobin and 
catalase. 
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preliminary, the separations observed in some systems are quite good 
(see hemoglobin in figures 15-17 and catalase in figure 17) , and 
indicate that it is possible to focus proteins on a pH gradient 
formed with only two or three simple ampholytes. 

4.2. Direction of Continuing Research. 

for simple ampholyte IEF to be of maximum, utility in protein 
separations, three types of focusing systems should be developed. In 
order of developmental priority they are: first, narrow-range focusing 

systems, which collectively span the pH range from 2 to 11 and indivi- 
dually cover a single pH unit utilizing a small number of simple 
ampholytes; second, a broad-range focusing system, which by itself 
covers the pH range from 2 to 11; third, a flexible method of designing 
custom-tailored pH gradients to fit specific and unusual separation 
problems. Our goal is to develop all three types of systems in forms 
suitable for both analytical and preparative application. 
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The objective of the present study is to develop a mathematical description of 
an Ampholina-free, multicomponent, isoelectric focusing system in the steady state, 
i.e., to obtain the means of computing ampholytes’ concentrations, the pH, conductance, 
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INTRODUCTION 

In isoelectric focusing (IEF) a stable pH gradient is created by 
electrophoretic migration of carrier ampholytes to stationary positions 
in the separation column. Ampholyfeic sample components introduced in the 
system distribute themselves on the pH gradient so that they exhibit con- 
centration maxima at their respective isoelectric points. In regards to the 
processes which determine their distribution on the column, the distinction 
between carrier ampholytes and sample ampholytes is purely semantic. 

1 

The existing theory of steady-state IEF is attributable to Svensson 
and was developed In conjunction with the introduction of the technique. 
Svensson addressed the one-component system, i.e., the concentration dis- 
tribution of a pure end component; and also the concentration distribution of 
an ampholytic component added to a system whose pH and conductance are known 
as a function of the distance along the separation axis. Such systems were 

rendered possible by Vesterberg's development of a random mixture of carrier 

2 3 

ampholytes , commercially introduced under the trade name of Ampholine. 

Their ready availability assured the remarkable grox^th of the various IEF 

techniques, but also rendered unnecessary further development of the theory 

of steady-state IEF. There have been no significant contributions to this 

theory since Svensson. 

Recent attempts to develop alternate means of forming stable natural pH 

4-7 

gradients have Been based on a purely empirical approach. We feel that 
this goal -would he furthered by the development of an explicit mathematical 
theory of multiple ampholyte IEF. This paper represents our approach in this 
direction, and we hope that it will contribute to an understanding of the 
theoretical Basis of IEF. 


MODEL FORMULATION 


The construction o£ a mathematical model of IEF is based upon the 
following concepts: 

a. the concentration of component subspecies are described by 
equations of chemical equilibria, 

b» in the steady state, a balance exists between the mass transports 
resulting from electromigration and from diffusion, 

c. the condition of electroneutrality prevails as a first approxi- 
mation, and 

d. the current density is constant throughout the separation column. 
The specific goal of the present theoretical treatment is to develop a mathema- 
tical description of a multicomponent IEF system at every point along its axis, 
i. e., to obtain the component concentrations, pH, conductance, and related 
parameters specified in terms of the components' mobilities, diffusion coeffi- 
cients, and proton dissociation constants. 

The theoretical treatment presented here assumes a good measure of ideality, 

“ft 

specific assumptions were as follows: (1) a negligible contribution of H 

and OH to the conductance within the pH range of interest, (2) only the two 
apparent dissociation constants CpK's) nearest the isoelectric point (pi) of 
an ampholyte exert significant influence upon its behavior in the steady state, 
(3) component subspecies have only a single absolute ion mobility, and (4) 
electroosmosis and radial and longitudinal temperature gradients are not 
significant. 

FUNDAMENTAL RELATIONS 

, The equations have been developed far both the two and three-component 
.systems. The three-component system is used to illustrate the fundamental 
relations. The concepts and treatment are generally applicable, and generalized 
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forms of the critical equations are given Below in the discussion of the 
method used in solution of the equations in computer simulation. 

The concentrations of component subspecies can Be expressed in terms of 
the total concentration of that component and the hydrogen ion concentra- 
tion h: 


h 2 C, 


'<* +' hK^ + K n K i2 ) 


Cl) 


c - = K ll K i2 C i 

1 (h 2 + hK^ + Vl2 ) 


( 2 ) 


^il c l 

a 2 + KK^ + Kil E 12 ) 


( 3 ) 


where i may Be a, b, or c representing the three components of the system; i 
is the concentration of cationic ampholyte; C^, the anionic ampholyte; and 
the uncharged and zwit ter ionic ampholyte. K_^~ and are the first 
dissociation constants of the ampholyte when titrated from the isoelectric 
state with acid and Base respectively. 

At any point the net mass transfer of each component resulting from 
( electromigration can Be described in terms of the charged subspecies concen 
trations as 



e 


I /-+ 

ir u i (c i 



( 4 ) 


where 1 is the current (amperes) ; 
conductance (mhos cm 2 ) ; and u^. 


2 

q, the cross-sectional area (cm ] ; ft, the 

2 -1 —1 

the component mobility (cm volt sec ) , 


H* + 


The conductance is given by 

K =* F/10G0 { u a (G^ + C") + u^CC* + CP' + u c (C+ 4* c“) } (5) 

where F is the Faraday constant and the component subspecies concentrations are 
molar. 

At any point the diffusional mass transport of each, component is described 
as 

J d “ * D i C6) 

2 —i 

where is the. diffusion coefficient Ccm sec 1 and x is the independent 
variable representing the distance along the column. Since in. the steady state 
the net mass transport of a chemical component resulting from electromigration 
must balance the. transport from diffusion, we can write 

- -jJ (7} 

e d 

and using equations 1 <3-1 and Q51» this becomes 


dc./dx => ~ — u. cct - c7) 

V D^qiC 1 i 2 - 


(3) 


Taking Into account equations Cl) and C2) we obtain explicit form of the equations 


dC /dx 

a 


dC. /dx 
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dC /dx 
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1000X 
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D qF 

(h 2 + »K al+ Kal K a2 3 


1000X 
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1000X 
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(9) 
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where the function $ has the form: 

u a C a (h 2 + K al K a2 ) + (h 2 + + u g C a (h 2 * K el K e2 > (9a) 

(h 2 + lyc ai r S a iS a2 ) Ch 2 + (h 2 + hS cl + s al K o2 ) 

The electroneutrality condition for the three-component system may be 
stated as 

h + C + + ct + C + = K /h + c“ + c7 + c" (10) 

a b c w a d c 

where is the ion product of water. 

NUMERICAL SOLUTION AND COMPUTER IMPLEMENTATION 
We present briefly the discretization procedure of the model equations, 
the numerical Integration technique, and the computer implementation of the model. 

the differential equations for C^, and 

W C b’ C c’ h > 

f b^ C a’ C h’ C c’ (11) 

f c tG a’ C V C c a h) 


It is convenient to write 
concentrations in the form: 


10001 

dC a /te * dIT 


dc /Sx = ,ioopi 
c b' < “ d.,f 


10001 

dC / ds ~ — — =■— 
c' D qF 


where f , f t , f are the functions specified in relation C3I. The functions 
a b c 

£ , f^, f c satisfy the Lipschitz condition with respect to all arguments. 
Using equations (X) and (2) we can write equation (10) in the form of a 


polynomial: 


8 

£ 

j.^o 



= 0 


( 12 ) 


where b g =* 1. 
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The coefficients b . of the polynomial are linear functions of C , C. , C 

2 3L D C 

of the form 


b* ■» b,_, + b. C +'b-. C-+b 4 C 
2 jF ja a jb b jc c 


( 13 ) 


and the parameters b^, b^ a , b^^, b^ depend on the dissociation constants 
K ah* ^k* K ch (^l*^) > and the ion product of water K^. These coefficients 
are presented explicitly in the Appendix* The solution of system (11) (12) , 

g 

which consists of differential and algebraic equations , is in general 

not unique* The uniqueness is guaranteed by picking one specific root of 

the polynomial. This has to be based on physical considerations. 

The differential equations (11) were discretized using the Runge-Kutta- 
9 

Merson scheme t which can be described as follows: 

C ,, =* C + Ax il) (G ,C. ,C .h ,Ax ) 
an+1 an n v an* bn J cn* V n' 


C bn+1 " C bn + ^n % ^an’Sn^cn’V^V 


(14) 


C _ . , =* C + Ax ijf (C ,C. ,C ,h ,Ax ) 
cnhl cn n r c an 3 bn’ cn J n J n' 


where 

Ax * x ,, - x is the integration stepsize. 
n nrl n. 

The functions ifj , , i{/ specified by the chosen discretization are computed 

SSL D C 

using the right sides of equations (11) . The subscript n added to the 

variables in the above aquation indicates the discrete values corresponding 

to the value of the independent variable x 

n-1 

x » x -f* F Ax, (15) 

a 0 1 

where x is the initial value, 
o 



The coefficients b . of the polynomial (12) are computed for each step x 

3 11 

based upon the values C , C, and C . Then the roots are extracted using 

the Newton-Raphson iterative technique. Having computed the roots we select 

the real positive root(s). When there is more then one positive real root, 

the root closest to the former value of h , is substituted for h . When 

n-1 n 

there is no real positive root the former value of h ^ is retained. This 
procedure can he described as follows: 

h “ min {r^J Re(r^) > 6 A Tm(r^) = 0} (16) 

Wi 

where Re(r^) is the real part of the root r^, and Im(r^) is the imaginary 
part of the root r^(Tsr=l,2, . . . ,8} ; h^ * ^ n _i the above set is empty. 

The program for numerical integration of system (11) (12) was prepared 
using FORTRAN based simulation software DAREP"^ available at the University 
of Arizona. The program was Implemented on the CYBER 175 digital computer. 
The polynomial (12) is solved by- the subroutine called at each integration 
step. The subroutine extracts the roots, tests them and substitutes appro- 
priate value for h^ following the formula (16) , When there are multiple 
real positive roots or there is none, an appropriate message is printed out. 
The integration is performed from the initial value x q until an end value 
x ™ x ma ^ . The step size Ax^ is automatically controlled by the Runge-Kutta- 

Merson integration subroutine to maintain local truncation errors within 

-3 -5 

assigned bounds. We assumed relative errors of 10 - 10 . In the cases 

studied the polynomial displayed only one positive real root which was 
therefore the only root realistically acceptable. The selection of this 
root of the polynomial at each integration step guarantees the uniqueness 
of the solution. 
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RESULTS' 


The model was exercised, at the level of the two-component system hy 
using values for mobilities, diffusion coefficients, and dissociation constants 
representative of glutamic acid and histidine* The constants assumed in the 
calculations are reported in Table X. The separation column was assumed to 
be 25 cm in length and of 0*79 cm cross-sectional area. The anode was 
positioned at the top of the column and the cathode at the bottom. The pre- 
dictions of the model for the two-component system are illustrated* in 
Figures 1, 2, and 3. 

The model was also exercised at the level of the three— component system 
by again using glutamic acid and histidine, and introducing a third component 
with an intermediate isoelectric point* The predictions of the model for 
the three-component system are illustrated in Figures 4, 5, and 6* 

DISCUSSION 

This work represents preliminary steps in the development of mathematical 
model for IEF. The results presented indicate that the method of solution 
chosen is capable of being used successfully. We are in the process of experi- 
mentally verifying the model’s predictions, and intend to perform parameter 
identification and remodeling as required. Although we do not believe the 
assumptions presently employed are unreasonable for preliminary theoretical 
consideration of the system, a more concrete position may not be assumed prior 
to experimental evaluation. 

The model is being generalised to consider an indefinite number ox 
components. At the level of the n-component system, it will provide a basic 
tool for the design and optimization of natural isoelectric focusing gradient 
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systems 


Recently a great deal a£ Interest has focused upon the modification 

4 5 

of pH gradients established by commercially available carrier ampholytes * 

and the generation of natural pH gradients by use of simple mixtures of 

6 7 

buffers * . It is generally appreciated that the quantities of individual 

ampholytes utilised and their specific characteristics must uniquely determine 

the shape of the pH gradient established, the sample capacity, and the 

resolution of separation. It should also he appreciated that the empirical 

approaches presently employed suffer from the absence of a mature theory 

of IEF. Significant progress cannot be reasonably expected without the 

development of a valid theory to accompany the empirical thrust. 

The computer simulation presented above has shown that the characteristics 

of gradients formed using a limited number of well-defined ampholytes will 

be far more dependent on their concentrations and the field applied than 

is the case when Ampholine is used. The reduction to practice of such novel 

buffer systems will be facilitated by the development of the automated data 

14 

collection apparatus described in another paper in this Volume . 

This research was supported in part by NASA, grant NSG-7333. 
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APPENDIX 

The coefficients of the polynomials are as follows: 
b l = - K a l K a 2 K t,l K b2 K cl K =2 K W 

K al K a2 K bX K b2 K cl K c2 CC a + ' C b + V 

CK al K a2 K bl K b2 K cl + K al K a2 K bl K cl K c2 + K al K bl K b2 K cl K 0 2) K » 

V- K al K a2 E bl K b2 K C A2 

CK al K a2 K bl K b2 K cl + K al K a2 K bl K cl K c2 } C a 
t ' K al K a2 K bl K b2 E al + K al K bl K b2 K cl K c2 ) C b 

- < K aAA2 K cl K c2 + K al K a2 K bl K c 1 K c2> C c 

CK al K a2 K bl K b2 + + K al K bl K b2 K cl + K al K a2 K cl K c2 

+ K al K bl K cl K c2 + K bl K b2 K cl K c2 ) K -j 

h ‘ Vzlhlhlh 1 + K aX K a2 K bl K cl K a2 + K aAlWc2 

+ CK 51 K b2 K cl K c2 * *«A*l*b2 - K al K a2 K cl K c2 * ‘u’tfuV C a 

+ <K al K a2 K al K c2 ~ K al K a2 K bl K b2 " K bl K b2 K cl K c2 " K al K bl K b2 K cl ) C b 

+ «a lVt*bl*b2 * K al K a2 K al K c2 " hlhlhlhl ' K aAl K cl K c2> C c 

- ^al^l + VA + E al K a2 K al + K aAl K cl + ^blVcl 

+ K al K cl K c2 + Wc2> % 

b 5 “ K al K a2 K bl K b2 + K AVd + E ilWd + K al K a2 K al K c2 

+ S alWc2 + K bl K b2 K cl K c2 

+ (K B1 K B2 K =1 + hlhlhz - K al K a2 K cl ' C a 

+ ®al K a2 K cl + K al K =l K c2 ' hlhl K cl ~ K al K bl K b2> C b 

+ + £ aAA2 " K al K ci K a2 “ K bl K cl K c2 ) C c 

+ «al K a2 + \lhl + + K al K cl + Vd + K =l E c 2> K w 


- 11 - 


+ + K al K a2 K cl + + ^1^2^01 

+ Wc2 + *BlWc2 


1 

1 

<%*2 

j. 

K =l K a2 

+ hihi 

- K al E a2> 

G a 

+ 

X 1 X 2 

+ 

K cl K c2 

+ K aAl 

" XlXP 

C b 

+ 

Xl K a2 

+ . 

Wl 

+ K aAl 

' K cl K c2> 

C 

C 

- 

Xl + % 

* K cl> 

K 

w 



B 7 ^ 

K al 5 a2 

+ XAi + hihz 

+ K ,K , • 

al cl 

H XAl + E el K c2 

+ 

Xl + 

Xl 5 C a + 

Xi + 

K al> C b + 

Xl + *bl> C = 


K. 







vr 






b 3 = 

K - + 

al 

^1 

+ K , 
cl 

4* G + 
a 

°b + °c 



- 12 - 


Table 1: Assumed Constants* 


Component 

P% 

PK 2 

5 2-1 -1 
u(10 )cm V sec 

D(10^)cm^sec 

A - Glutamic acid 

2.19 

4.25 

31.7 

0.80 

B - Histidine 

6.00 

9.17 

32.4 

0.73 

C — Unspecified component 

4.00 

6.00 

30.0 

0.80 


* pKL Values were taken from reference 11; diffusion coefficient for histidine 

from reference 12; diffusion coefficient for glutamic acid was calculated 

13 

from the values of the equivalent conductance ; mobilities were taken as 
an average of the values for the anionic and cationic species. 
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Fig. 1: The concentration profiles of the two-component 

IEF system using glutamic acid and histidine. The amperage 
employed was 4x1 0 -7 A. 
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Fig, 2: The pK profile of the two-component system. 

The conditions were identical to Fig. 1. 
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Fig. 3: TTxe conductance profile of fcfie two-component 

system. The conditions were identical to Fig. 1. 
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system. The conditions were identical to Fig. 4 
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Fig. 6: The conductance profile of the tBree- 

component system. The conditions were identic?,! to 
Fig. 4. 
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MATHEMATICAL MODEL OE STEADY-STATE ISOELECTRIC FOCUSING 
SYSTEM OF N-COMPONENTS 


G. A, P ALUS INSKI 
Title 

Mathematical Model of a Steady-State Isoelectric Focusing 
System of N-Components 

Abstract 

The paper presents basic equations expressing the 
concentrations of cationic and anionic ampholytes in terms of 
hydrogen ion concentration, corresponding dissociation con- 
stants and ampholyte concentrations . Balance between electro- 
migration and diffusion movement determines the set of 
differential equations. Solution of these equations gives the 
ampholyte concentration as function of distance x, along the 
column axis. The differential equations are coupled by the 
polynomial of order 2N+2, which results from the electro- 
neutrality equation. The coefficients of the polynomial are 
complicated functions of dissociation constants and ampholtye 
concentrations. The recurrence formula for these coefficients 
is derived in the first part of the paper. 
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Introduction 

In an earlier paper ( ) a mathematical model was 

presented for a three-compcpnent isoelectric focusing system. 
The model provided the means of computing ampholytes' con- 
eentrations , the pH and the conductance at every point along 
the column axis in terms of the components' mobilities, 
diffusion coefficients, and proton dissociation constants. 

The model consists of an ordinary differential equation for 
each component and a single algebraic coupling equation. The 
differential equations express the steady-state balance that 
exists between the mass transports resulting from diffusion 
and electromigration. The algebraic polynomial is a statement 
of the condition of electroneutality . The equations of the 
two and three component models were derived directly from the 
fundamental relations . In order to conveniently implement the 
model with higher multiple component systems, the model 
equations have been generalized and an algorithm has been 
developed to generate the differential equations and coupling 
polynomial for an N component system. 
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Notation 

_L, 

C. ~ the concentration of cationic ampholyte (ith comp.) 


Ct - the concentration of anionic (units 
h - hydrogen ion concentration (units) 


K il' K i2 


the first dissociation constants of the ampholyte 


(titrated from isoel. state with acid and base 
respectively) . 
u^ - component mobility 

q - cross-sectional area 

-- diffusion coefficient 
K - ion product of water 


-B3- 


1 


1. BASIC EQUATIONS 


+ h2 "i 

C. = i 


h + hK. 1 + K. 1 K i2 


C i ,2 


K, i K. „ C. 
xl x2 x 


h + hK. , ■ + K. -K. n 
xl xl x2 


1 = 1/ 2, . . . ,N 


CONDUCTANCE 


N 


K = lfoo J 1 u i (C i + + C i 


ELECTROMIGRATION 


je - I? u i (c i + - c i"> 


DIFFUSION 


dC. 


Jd -- D. 


i dx 


x x 

BALANCE (Je = -Jd ) 


ac^ 

dx 


Dj^qK U i (C i " C i 3 


ELECTRONEUTRALITY 


N 

h + 2 C* 

i=l x 


+ N 

: K /h + 2 C. 

w 1=1 1 
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( 1 ) 



(3) 


(4) 


(5) 



2 


2 . DIFFERENTIAL EQUATIONS 


Cl)/ (2) + (3) 


N h + K.-K. 

K = F y T , xl x 

X000 "i h 2 K K . 

xl xl x2 


(1)/ (2) -> (4) 


dC. _ h 2 - K../K.- 

x _ I u xl x2 

dx ^qK x h 2 + hK.' + K.-K. , 

xl xl x2 


(6). + (7) 


dC. 10001 u. h - K. -K. 0 -t 

X __ X xl x 2 -L Q 

ax “ FD i q h2 + “ii + K ii K i 2 5 1 


N h‘ + K ± K 

Q = S ^ C. 

i“l h + “ll + K il K i2 


3 . POLYNOMIAL 


(1)/ (2) + (5) 


» H 2 =i „ »' K il K i2 C i 

h + S -5 = = K 7h + S - s - — — 

i=l h 2 + hK u + K. lKi2 W i=l h + hK^ + K ilK . 2 


h 2 + l h C l~ K il K 12 hC i 

i=l h 2 + hK. x + K i;L K i2 


= ^ + hK ji + K ji K j2> (11) 
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3 


h 2 ! (h 2 + h Kjl+ K. 1 K j2 ) + I (h 3 C i -hK. 1 K. 2 C i ) ? (h^hK.^K.^ 


j=l 


j=l 

j+i 


* 2 

= K ¥ (h + hK., + K . ,K . „) 


w 


3=1 


ji *'jr ]2' 


( 12 ) 


Notation: 


A. = K. 
i ll 


B. = K.-K.„ 

i il i2 


(13) 


~ N - N - N P 

h ¥ (h +hA.+B.) + E (h C.-hB.C.) tt (h + hA. + b.) “ 
j-i 3 3 i=l 1 1 1 j=l 3 3 


• 1 * 
Jt 1 


N 

= K ¥ (h + hA. + B.) 
W j=1 3 3 


(14} 
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W 9 2 9 2 

P w = TT (h +hA.+B.) a (h +ha 1 -rB 1 ) (h +hA-+B-) (h‘ +hA-+B 0 ) 

JM ^ j J -LX /• si JJ 


. - (h + hA + B_ ) 
n n 


N 


1 " it (h + hA. + B.) 
1 j=2 ^ ^ 


N 2 

P N 2 = IT (h + hA. + B.) 
j=l ^ 3 

jf 2 


N 2 

" hA j +B j ) 

, j+i 


N 2 

P w w i = IT (h + hA. + B.) 
M, N-X 3 3 

j^N-1 


P H,H = (h2 + + B j> 


Using the notation (15) we can write (14) in the form 


(15) 


h2 + ji <1,3c i ' hB i C i’ Vi = K « P N 


(16) 


The (16) is the general form of the polynomial in the case of 
N- component. 



5 


4. THE RECURRENCE RELATIONS FOR THE POLYNOMIAL COEFFICIENTS 

4 . 1 Fundamental relation 

The polynomial (16) c;«.n be written in the form 


fa2 + 


» 3 

E (h C.-hB.C.) 
, r i r 


p 

Nr 1 




0 


(17) 


We shall use the notation 


V W “ + ' hB i C i> Vi - K w P N (1S ’ 

Expressing F and P , explicitly in (18 ) , performing 

JM N t — 

suitable operation .and grouping the terms with respect to the 
powers of h we obtain 


2N+2 . N 

* N (h) = Z h 1 a i (19) 

' s i=o 

N 

where . a i - coefficient of the polynomial in case of N-components . 
The polynomial for the case of N-fl components has the form 


h 2 P 


N+l 


N+l 


i=l 


(h 3 C ± 


- h B. 


C i> 


N+lri 


- K 


w N+l 


= 0 


( 20 ) 


derived directly from (17) . 

The formula (20) can be written in the form 


h 2 P 


N+l 


N 
+ Z 
i=l 


(h 3 C . -hB C. ) P__. 1 . -K P _ + (h 3 C *“hB - C._ , - ) 

1 2. N+l r l w N+l N+l N+l N+l 


P =0 

N+l r N+l 


( 21 ) 


-B8- 
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From (15) we obtain 


P N+1 P N' h * hA N+l + B N+1 J 


P,,,-, . — . (h 4- hAv. , + B T - ) 

N+l, x N,x N+l N+l 


i ~ 1 z 2 j . « * * ^ N 


? N+l, N+l P N 


Using (22) in (21) we obtain 


N 


h P N + . S . C^C.-hB.C.JP - Vn (h + “W + B N + 1> 
1=1 


+ 


+ (h C N+1 “ h E N+1 C Ntl ) ° 


Taking into account notation (23) we obtain the polynomial 
for N+l components expressed in terms of polynomial for N 
components 


^ ^ +hA N+l + B N+1 J + C N+1 ” hB N+l G N+l^ P N ° 


The above polynomial for N+l components can be divided into 
two parts 

2 


F* ^ + k ^KT-J.1 + - -NT-l-l ^ 


N+l N+l’ 


■S2 - (h C„.,. - h B . , c -p. 

N+l N+l N+l N 


-B9- 


( 22 ) 


(23) 

(24) 

(25) 

(26) 
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Substituting (19) to (25) we obtain after some manipulations 


2N+4 ■ 

I - 2 11 ( a 1 ._ 9 + 

i=o 


i-2 a i-l N+l i W+1 J 


with 


a 2N+2 1 


a -k a 2N+2+k ° 


for k = 1 , 2 ,... . 

In order to determine the coefficients of £1 we have to find 
the expression for 


4.2 Finding the recurrance relations for 


We shall use the notation 


2N . „ 

V = l h? N a. 
H j=o ^ 


) a 2N 


(28) 


From (22) we have 


P N+1 P N (h + h A N+1 + W 


(29) 


Substituting (2£„ to (29) we obtain after some manipulations 


2N+2. • XSI N 

P N+1 " . 2 h ( a j-2 + a j-l A N+1 + a j B N+1 } 

3-0 


( 30 ) 


where 


a -k Na 2N+k 0 


for k = 1/2; 
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The formula (30) permits computing the coefficients of 
based on the coefficients of 


N+l 

ct . 


N , N . .. 55 _ ' 

a* 0 + a . - A . ■+ a. B _ 
j— 2 j-1 N+l n N+l 


this way we can successively compute P^/- * * * ■/ etc 


Substituting (28) to (26) yields after some manipulations 


_ 2 ^+4 • N N 

" , S o h ( a i-3 “ a i-l E N+1 } C N+1 
1=0 


(32) 


N a_ k = N ct 2N+k = 0 for k = 1,2,3 


N 


a 2N 1 


4«3 Recurrence notation for the polynomial coefficients . 

The polynomial for N+l components given by formula (24) 
can be written in the form 


w (h) = r+ n < 33 > 

where f" and S2 are specified by (25) , (26) or (27) , (32) . 
Substituting (27) and (32) to (33) yields 




N+l 


(h) = 


2N+4 
2 3 
i=o 


i N 

C 


a i-2 


N 


i-1 N+l 


, N _ , ,N N _ , p 

+ a i B N+l + a i-3“ a i~l B N+l C N+1 


3 


(34) 


-Bll— 


From (34) we obtain the recurrence formula for the 
coefficients of the polynomial in the form 


N+l _ N , N , ,N N . _ 

a l a i-2 + a i-l ^N+l + N a i B N+l + ^ a i-3 a i-l B N+l^ C N- 


i = 0/1,2,..,, 2N+4 


:h+i 

{35} 


N _ N n 

-k a 2N+2+ k U 


N = N 
a -k a 2N+k 


= 0 


k - 1,2,3 


-B12- 


APPENDIX C 


INDEX OF LIBRARY OF COMPUTER SIMULATIONS 



Identification 

Dissociation 

Constants 

C-l 

C-2 

Crt, 

Concentrations 

C 

Hobilities 

c. 

Diffusion Constant 

Ho . Date 

Tima 

A1 

A-2 

B-l 

8-2 

A 

D 

A 

8 

A 

n 

C 

1. 4/7 

14 : 37 

2.19 

4.25 

6.0 

9.17 



1«- 5 

3.4x-3 

3.4x-3 


30. 0x-5 

30. Ox-5 


0.9x-5 

0.9X-5 


2. 4/7 

14:03 

2.19 

4.25 

0.0 

9.17 



lx- 5 

2.78-6 

2.70X-5 


30. Ox-5 

30. 0x-5 


0.9x-5 

0.9X-5 


3 . 4/12 

9:39 

2.19 

4.25 

6.0 

9.17 



lx- 6 

7.2x-3 

8.1x-5 


31.7x-5 

32.4x-5 


0.8x-5 

0.73X-5 


4. 4/12 

9:48 

2. 19 

4 , 25 

6.0 

9.17 



lx-7 

7.2x-3 

8.1x-5 


31.7x-5 

32.4x-5 


O.Ox-5 

0.73x-5 


5* 4/12 

10:10 

2.19 

4.25 

6.0 

9.17 



4x-7 

7.2X-3 

8.1x-5 


31.7X-5 

32 ,4x-5 


O.Cx-5 

0.73X-5 


6. 4/12 

12: 21 

2.19 

4.25 

6.0 

9.17 



2x-7 

7.2z-3 

B.lx-5 


31.7X-5 

32.4X-5 


O.Sx-5 

0.8x-5 


7. 4/12 

9:53 

2.19 

4.25 

6.0 

9.17 



5x-7 

7.24X-3 

0.1x-5 


31.7X-5 

32.4X-5 


0.8X-5 

0 . 73x-5 


0. 4/12 

12:33 

2.19 

4.25 

6.10 

9.17 



2x-7 

5.46x-3 

1.23X-3 


31.7X-5 

32.4x-5 


O.Bx-5 

0.73X-5 


9. 4/12 

12:39 

2.19 

4,25 

6.0 

9.17 



5x-7 

7.69X-8 

7.27x-3 


31.7X-5 

32,4x-5 


O.0X-5 

0.73x-5 


10. 4/12 

12:43 

2.19 

4.25 

6.0 

9.17 



4x-7 

7.24X-3 

Q.lx-5 


31.7X-5 

32.4X-5 


O.Ox-5 

0.73x-6 


11. 4/12 

12:26 

2.19 

4.25 

6.0 

9.17 



4X-7 

7,24x-3 

8.1x-5 


31.7X-5 

32.4x.-5 


O.Bx-5 

0.73X-5 


12. 4/24 

21:23 

2.19 

4.25 

6.0 

9.17 



4x-7 

7.24x-3 

O.lx-5 


31.7X-5 

32.4x-5 


O.Bx-5 

0.73x-!> 


13. 4/2G 

16:44 

2.19 

4.25 

6.0 

9.17 



4x-7 

7.24X-3 

a.lx-5 


31.7x-5 

32.4X-5 


O.Dx-5 

0,73x-4 


14. 5/2 

16:16 

0.19 

6.25 

4.11.17 



4x-7 

7.24X-3 

a.ix-s 


31.7X-5 

32.4x-5 


O.Ox-5 

0.73x-5 


15. 5/2 

17:00 

0.19 

6.25 

4.0 

11:17 



4x-7 

7.24x-3 

8 , lx- 5 


31.7X-5 

32.4X-S 


O.flx-5 

0.73X-5 


16. 4/29 

21:16 

2.19 

4.25 

6.0 

9.17 

4.0 

0.0 

4x-6 

7.24X-3 

B.lx-5 

6. Ox-4 

31.7X-5 

32.4X-5 

30. 0x-5 

O.Bx-5 

0.73X-5 

O.Ux-5 

17. 4/30 

12:16 

2,19 

4.25 

6.0 

9.17 

4.0 

6.0 

4x-6 

7.24x-3 

B.lx-5 

6. Ox-4 

31.7X-5 

32.4X-5 

30. 0x-5 

O.Bx-5 

0.73X-5 

O.Ox-5 

10. 4/30 

12:18 

2.19 

4.25 

6.0 

9.17 

4.0 

6.0 

4x-6 

7.24X-3 

B.lx-5 

6. Ox-4 

31..7X-3 

32.4X-5 

30. Ox-5 

O.Bx-5 

0.73X-5 

O.Bx-5 

19. 4/30 

12:37 

2.19 

4.25 

6.0 

9.17 

4.0 

6.0 

4x-7 

7.24x-3 

B.lx-5 

6. Ox- 4 

31.7x-5 

32.4x-5 

30. Ox-5 

O.Ox-5 

0.73X-5 

O.flx-5 

20. 5/2 

16:40 

2.19 

4.25 

6.0 

9.17 

4.0 

6.0 

6x-7 

7.24x-3 

fl.lx-5 

6. Ox- 4 

31.7x-5 

32.4x-5 

30. Ox-5 

O.Bx-5 

0.73x-5 

O.Bx-5 

21. 5/2 

10:55 

2.19 

4.25 

6.0 

9.17 

4.0 

6.0 

5x-7 

?.24x-3 

B.2x-5 

6. Ox-4 

31.7X-5 

32.4X-5 

30. Ox-5 

O.Bx-5 

0.73X-5 

Q.BX-5 


*% 

o v 




Identification 

Dissociation 

Constants 

C-l 

C-2 

Crt. 

Concentrations 

Mobilities 

Diffusion Constant 

Mo . Date 

Time 

A1 

A-2 

D-l 

B-2 

A 

B C 

A 

13 C 

A 

n c 

1. 6/15 

13:22 

2.0 

4.0 

6.0 

0.0 



4x-7 

7.24X-3 

O.lx-5 

30. Ox-5 

30. Ox-5 

O.Bx-5 

O.Ox-5 . 

2. 6/15 

13:23 

2.0 

4.0 

6.0 

8.0 



4x-7 

7 . 24x~3 

8. lx- 5 

30 , Ox-5 

30. Ox-5 

0.0x-5 

O.Ox-5 

3. 6/15 

13:24 

2.0 

4.0 

5.0 

9,0 



4x-7 

7.24X-3 

B.lx-5 

30 . Ox-5 

30. Ox-5 

O.Dx-5 

O.Ox-5 

4. 6/15 

13:28 

2.0 

4.0 

5.0 

9.0 



i!Jc-7 

7.24X-3 

B.lx-5 

30. Ox- 5 

30. Ox-5 

Il.llx-G 

n.ox-b 

5. 6/15 

13:42 

1.0 

5.0 

6.0 

8.0 



4x-7 

7.24X-3 

B.lx-0 

30. Ox-5 

30. Ox-5 

n.flx-5 

O.Bx-5 

6. 6/15 

13:32 

2.0 

4.0 

4.0 

10.0 



4x-7 

7.24x-3 

B.lx-5 

30. Ox-5 

30. Ox-5 

0,3x-6 

O.Ox-5 

7. C/15 

13:43 

3.0 

5.0 

6.0 

8.0 



4x-7 

7.24X-3 

B.lx-5 

30, Ox-5 

30. Ox-5 

G.Bx-5 

O.Bx-5 

0. 6/15 

13:45 

4.0 

6.0 

6,0 

8.0 



4x-7 

7.24X-3 

B.lx-5 

30. Ox- 5 

30, Ox-5 

O.Bx-5 

0. Ux-b 

9. 6/16 

9:47 

5.0 

7.0 

6.0 

0.0 



4x-7 

7.24X-3 

8.1X-5 

30. Ox- 5 

30. Ox- 5 

9, Ox- 5 

O.Ox-5 ■ 

10. 6/16 

9:50 

5.5 

7.5 

6.0 

8.0 



4x-7 

7.24X-3 

B.lx-5 

30. Ox-5 

30. Ox-5 

O.Bx-5 

n . Bx- 5 

U, 6/16 

9:51 

2.0 

4.0 

6.0 

0.0 



4x-7 

7.24X-3 

8. lx- 5 

30. Ox-5 

30. Ox- 5 

O.Bx-5 

n.Flx-5 

12. 6/16 

9:59 

1.0 

9.0 

6.0 

8.0 



4x-7 

3. 0x-3 

3. Ox-3 

30. Ox-5 

30. Ox-5 

O.Ox-5 

O.Ux-5 

13. 0/16 

10:00 

2.0 

8.0 

6.0 

8.0 



4x-7 

3. Ox-3 

3. Ox-3 

30. Ox-5 

30.0x.-5 

O.Bx-5 

n.Bx-5 

14. 6/16 

9:57 

3.0 

7.0 

6.0 

8.0 



4x-7 

3 , Ox- 3 

3. Ox- 3 

30. Ox-5 

30.0X-5 

O.Dx-5 

O.Rx-5 

15, 6/16 

9:54 

4.0 

6.0 

6.0 

8.0 



4x-7 

3. Ox-3 

3. Ox-3 

30, Ox-5 

30. Ox- 5 

O.Bx-5 

O.Bx-5 

16. 6/16 

9:53 

2.0 

4.0 

6.0 

8.0 



4x-7 

3. Ox-3 

3 . Ox- 3 

30. Ox-5 

30. Ox-5 

O.Bx-5 

o,nx-5 

17. 5/30 

3:49 

5.5 

7.5 

6,0 

8.0 



2x-6 

3. Ox— 3 

3. Ox-3 

30. Ox-5 

30. Ox- 5 

O.Bx-5 

n.0x-5 

IS. 6/30 

3:44 

1.0 

9.0 

6.0 

0.0 



4x-B 

1.5x-2 

l.Sx-2 

30. Ox-5 

30, Ox- 5 

0. Ox-5 

O.Ox-5 

19. 6/29 16:55 

5.0 

7.0 

6,0 

8.0 



4x-6 

3. Ox-3 

3. Ox-3 

30. Ox-5 

30. Ox-5 

O.Ox-5 

O.Bx-5 

20. 6/30 

0:52 

5.5 

7.5 

6.0 

8.0 



4x-6 

3. Ox-3 

3. Ox-3 

30. Ox-5 

30. Ox-5 

O.Sx-5 

O.Ox-5 

21. 6/23 15:43 

1.0 

9.0 

6,0 

8.0 



4x-6 

3. Ox-3 

3. Ox-3 

30. Ox-5 

30. Ox-5 

O.Bx-5 

O.Ox-5 

22. 6/30 

8:25 

1.0 

9.0 

6.0 

8.0 



4x-6 

3. Ox- 2 

3. Ox- 2 

30. Ox-5 

30. Ox-5 

O.Bx-5 

O.Ox-5 

23. 6/29 16:53 

5.0 

7.0 

6.0 

8.0 



2x-6 

3 , Ox-3 

3. Ox-3 

30. Ox-5 

30. Ox-5 

0.8x-5 

O.Bx-5 

24. 6/30 

3:22 

1.0 

9.0 

6.0 

8,0 



4x-6 

3. 0x-2 

3. Ox- 2 

30, Ox-5 

30. Ox-5 

O.Bx-5 

0 . Bx- 5 

25. 6/29 17:06 

1.0 

9.0 

6.0 

8.0 



4x-7 

3. Ox-3 

3. Ox-3 

30. Ox-5 

30. Ox-5 

O.Bx-5 

O.Bx-5 


ORIGINAL! PAGE IS 
OF POOR QUALITY 



Identification 

Dissociation 

Constants 

C-l 

C-2 

Crt. 

Concentrations 

C 

Mobilities 

C 

Diffusion Constant 

Ho. Date 

Time 

AX 

A-2 

0-1 

B-2 

A 

0 

A 

B 

A 

B 

C 

26, 6/29 

17:00 

2.0 

8.0 

6.0 

0.0 



4x-7 

3. Ox-3 

3. Ox-3 


30 , Ox-5 

30. Ox-5 


Q.Bx-5 

0. Ox- 5 


27. 6/29 

16:58 

5.5 

7.5 

6.0 

8,0 



4x-6 

3. Ox-3 

3. Ox-3 


30. Ox-5 

30, Ox-5 


O.Ox-5 

O.Bx-5 


28. 6/29 

17:02 

5.5 

7.5 

6.0 

0.0 



2x-6 

3. Ox-3 

3 . Ox-3 


30. Ox-5 

30. Ox-5 


O.Qx-5 

0, Ox- 5 


29. 6/30 

8:07 

5.0 

7,0 

6,0 

0,0 



2x-6 

3. Ox-3 

3. Ox-3 


30. Ox-5 

30. Ox-5 


O.Qx-5 

O.flx-5 


30. 6/30 

0:05 

5.0 

7.0 

6.0 

0,0 



<ix-6 

3. Ox-3 

3. Ox-3 


30. Ox-5 

30. Ox-5 


O.flx-5 

O.flx-5 


31. 6/29 

16:51 

1.0 

9.0 

6,0 

0.0 



Ox- 6 

3. Ox-2 

3 . Ox-2 


30. Ox-5 

30. Ox-5 


O.Ox-5 

O.flx-5 


32. 6/29 

16:09 

1.0 

9.0 

6.0 

0.0 



4x-5 

1. 5x-2 

1.5x-2 


30. Ox-5 

30, Ox- 5 


O.Ox-5 

O.flx-5 


33. 6/16 

10:19 

1.0 

9.0 

6.0 

Q.O 



<lx-7 

3. Ox-3 

3. Ox- 3 


30. Ox-5 

30. Ox-5 


O.Ox-5 

O.Ox-5 


34. 6/29 

15:01 

1.0 

9.0 

6.0 

0.0 



2x-6 

3. Ox-3 

3, Ox- 3 


30. Ox-5 

30 . Ox-5 


O.Ox-5 

O.flx-5 


35. 6/29 

15; i)0 

2.0 

0.0 

6.0 

0.0 



<lx-6 

3. Ox-2 

3. Ox-2 


30. Ox-5 

30. Ox-5 


O.Ox-5 

O.Ox-5 


36. 6/29 

15:06 

2,0 

0.0 

6.0 

0.0 



<lx-6 

1.5X-2 

l.Sx-2 


30. Ox-5 

30. Ox-5 


O.Ox-5 

O.Ox-5 


37. 6/29 

15:30 

2.0 

0.0 

6.0 

0.0 



<ix-6 

3. Ox-3 

3. Ox-3 


30. Ox- 5 

30. Ox-5 


O.Bx-5 

O.flx-5 


30. 6/29 

15:35 

2.0 

0.0 

6.0 

0.0 



2x~6 

3. Ox-3 

3. Ox-3 


30. Ox-5 

30. Ox-5 


O.Ox-5 

O.flx-5 


39. 7/11 

10:57 

5.0 

7.0 

3.0 

11.0 

7.0 

9.0 

lx- 6 

3. Ox- 3 

9. Ox-3 

3. Ox-3 

31.7X-5 

32.4x-5 

30. Ox-5 

O.Bx-5 

0.73X-S 

O.flx-5 

<10. 7/11 

11:05 

5.0 

7.0 

3.0 

11.0 

7.0 

9.0 

2x-7 

3. Ox-3 

9. fly- 3 

3. Ox- 3 

31.7X-5 

32. Ox-5 

30. Ox-5 

O.Qx-5 

0.73X-5 

0 . Ox-5 

01. 7/11 

10: 50 

5.0 

7.0 

3.0 

11.0 

7.0 

3.0 

lx- 6 

3. Ox-3 

6. Ox-3 

3. Ox-3 

31.7X-5 

32. Ox-5 

30. Ox-5 

O.Bx-5 

0.73x-5 

O.Ox-5 

<i2. 7/11 

11:07 

5.0 

7.0 

3.0 

11.0 

7.0 

9.0 

2x-7 

3. Ox-3 

6. Ox- 3 

3, Ox-3 

31.7X-5 

32.<lx-5 

30. Ox- 5 

O.Bx-5 

0.73X-5 

O.Ox-5' 

<13- 7/11 

10:00 

5.0 

7,0 

6.0 

0.0 

7.0 

9.0 

lx- 6 

3. Ox-3 

12. Ox-3 

3. Ox-3 

31.7x-5 

3 2. Ox-5 

30, Ox- 5 

O.Ox-5 

0.73X-5 

0 . Ox -5 / 

j i 

00. 7/11 

10:02 

5.0 

7.0 

6.0 

8.0 

7.0 

9.0 

lx- 6 

3. Ox-3 

9. Ox- 3 

3. Ox -3 

31.7X-5 

32. Ox-5 

30. Ox-5 

O.Ox-5 

0.73X-5 

O.flx-jS ! 

U5. 7/11 

10:31 

5.0 

7.0 

6,0 

0.0 

7,0 

9.0 

lx-6 

3. Ox- 3 

6. Ox-3 

3. Ox-3 

31.7x-5 

32.HX-5 

30. Ox-5 

O.flx-5 

0.73X-5 

O.Ox-5 ! 

06. 7/11 

9:27 

2,0 

0.0 

3.0 

9.0 



<ix-5 

6, Ox- 3 

6. Ox-3 


30. Ox-5 

30, Ox-5 


O.Ox-5 

O.Ox-5 

\ F 

j 

07. 7/11 

9:25 

2.0 

0.0 

3,0 

9.0 



<)X-6 

6 , Ox- 3 

6. Ox-3 


30. Ox-5 

30. Ox-5 


O.Ox-5 

O.Ox-5 


00. 7/11 

9:07 

2.0 

8.0 

<1.0 

10,0 



<lx-5 

6. Ox-3 

6. Ox- 3 


30. Ox- 5 

30. Ox-5 


O.Ox-5 

O.Ox-5 

t 

09. 7/11 


2.0 

0.0 

n.o 

10.0 



Ox-6 

6. Ox- 3 

6. Ox-3 


30. Ox-5 

30. Ox-5 


O.Ox-5 

O.Ox-5 


50. 7/11 

9:13 

2.0 

0.0 

<1.0 

10.0 



Ox-G 

6. Ox-3 

6. Ox-3 


30. Ox- 5 

30. Ox-5 


0,8x-5 

O.Ox-5 




Identification 

Dissociation 

Constants 

e-i 

C“2 

Crt. 

Ho. Date Time 

A1 

A-2 

B-l 

B-2 

51. 7/11 10:05 

2.0 

0.0 

4.0 

10.0 



]x-6 

52. 7/11 10:20 

2.0 

0.0 

4.0 

10.0 



<lx~7 

53. 7/ii 9:05 

1.0 

3,0 

2.0 

4.0 



4x-7 

5*1 . 7/11 9:03 

5.0 

7.0 

6.0 

B.O 



4x-7 

55. 7/12 15: <H 

5.0 

7.0 

6.0 

8.0 

7.0 

9.0 

lx-6 

56. 7/12 15:46 

5.0 

7.0 

6,4 

6.6 

6.0 

8.0 

4h-7 

57. 7/13 11: <13 

5.0 

7.0 

6.0 

0.0 

7.0 

9.0 

lx- 6 


Cl 

-ft- 

1 


Concentrations 


Mobilities Diffusion Constant 


A 

D 

C 

A 

D 

C 

A 

B 

C 

6. Ox-3 

6, Ox-3 


30. 0x-5 

30. Ox-5 


O.Bx-S 

O.fix-5 


6. Ox-3 

6, Ox-3 


30. 0x-5 

30, Ox-5 


O.Bx-5 

0. fix-5 


3. Ox-3 

3. 0x-3 


30. Ox- 5 

30. Ox-5 


0 . fix-5 

O.Ox-5 


3. Ox-3 • 

3. Ox-3 


30. 0x-5 

30. 0x-5 


0.8x-5 

O.flx-5 


2. 0x-6 

3. Ox-3 

2. Ox-6 

30. 0x-5 

30. Ox-5 

30. Ox-5 

0 . Ox-5 

O.Ox-5 

O.Ux-5 

3,0x-3 

4. Ox- 3 

3. Ox-3 

30. Ox-5 

30. 0x-5 

30. 0x-5 

O.Elx-5 

O.fix-5 

O.Cx-5 

3.0x-3 

2. Ox-7 

4. Ox- 10 

30. 0x-5 

30.0X-S 

30. 0x-5 

O.Ox-D 

0.8X-5 

O.fix-5 




r • 


Identification 

Dissociation 

Constants 

C-l 

C-2 

Crt. 

Concentrations 

C 

Mobilities 

C 

Diffusion Constant 

Ho, Date 

Time 

01 

A- 2 

B-l 

B-2 

A 

B 

A 

D 

A 

B 

C 

1. 8/2 

21:40 

2.19 

4.25 

4.0 

6.0 

6.0 

9.17 

Gx-7 

8. 0x-3 

8 . Ox- 5 

1. Ox-10 

31.7X-5 

30. Ox-5 

32.4x-5 

0. Bx-5 

O.Bx-5 

0.73X-5 

2, 0/3 

21:20 

2.19 

4.25 

4.0 

6.0 

6-.0 

9.17 

2x-G 

8. 0x-3 

1. Ox-4 

1. Ox-14 

31.7H-5 

30. Ox-5 

32.4X-5 

0. Bx-5 

0, Bx-5 

0.73X-5 

3. 0/!| 

11: 55 

4.0 

6.0 

6,0 

0.0 

8.0 

10.0 

2.2x-6 

8 . Ox- 3 

1, Ox-4 

l.Qx-26 

31.7x- 5 

30. Ox-5 

32.4X-5 

O.0X-5 

O.Ox-5 

0.73X-5 

4. 0/3 

9:07 

2.19 

4.25 

4.0 

6.0 

6.0 

9.17 

lx- 5 

G.Ox-3 

0i Ox-5 

l.OK- 11 

31.7x-5 

30. Ox-5 

32.4x-5 

O.Ox-5 

O.Bx-5 

0.73X-5 

5, Ho Date 

2.19 

4.25 

4.0 

6.0 

6.0 

9.17 

lx- 6 

0.OX-3 

l.X-4 

l.Ox-11 

31.7x-5 

30. Ox-5 

32.4X-5 

O.Qx-5 

0.8x-5 

0.73X-5 

6. 0/4 

9:24 

2.19 

4.25 

4.0 

6.0 

6.0 

9.17 

2x-6 

8. Ox-3 

1, Ox-4 

1 ,0x- 18 

31.7X-5 

30. Ox-5 

32.4x-5 

O.0X-5 

O.Ox-5 

0.73X-5 

7. 0/4 

9:30 

2.19 

4.25 

4,0 

6.0 

6.0 

9.17 

2x-6 

8. Ox- 3 

1. Ox- 4 

1. Ox-26 

31.7X-5 

30. Ox-5 

32.4X-5 

O.0X-5 

O.Bx-5 

0.73X-5 

8. 0/4 

9:43 

2.19 

4.25 

4.0 

6.0 

6.0 

9.17 

1.6x-6 

B.Ox-3 

l.Ox-n' 

1. Ox-26 

31.7X-5 

30. Ox-5 

32.4X-5 

O.0X-5 

0,8x-5 

0,73x-5 

9. S/if 

10:31 

2.19 

4.25 

4.0 

6.0 

6.0 

9.17 

1.7x-6 

0. Ox-3 

1.0x-4 

] .Ox-26 

31.7X-5 

30. Ox-5 

32.4X-5 

O.Bx-5 

0,8x-5 

0,73x-5 

10. 0/4 

10:40 

2.0 

4.0 

4.0 

6.0 

6.0 

0.0 

1.7x-6 

8. Ox-3 

1. Ox-4 

l.Ox-26 

31.7X-S 

30 . Ox-5 

32.4X-5 

O.Bx-5 

0.8x-5 _ 

0.73X-5 

11. 0/4 

11:06 

2.0 

4.0 

4.0 

6.0 

6.0 

D.O 

1.7x-6 

0 . Ox- 3 

1, Ox-4 

1. Ox-26 

30. Ox- 5 

30 . Ox-5 

30. Ox- 5 

o’. Bx-5 

O.Gx-5 

O.0X-5 

12. 8/4 

11:20 

4,0 

6.0 

6.0 

0.0 

8.0 

10,0 

1.7x-6 

0 . Ox- 3 

1. Ox- 4 

1. Ox-26 

31.7X-5 

30 . Ox-5 

32.4x-5 

0.8x-5 

O.Gx-5 

0.73X-5 

13. 0/4 

11:34 

4.0 

6.0 

6,0 

0.0 

0.0 

10.0 

1, 0X- 6 

O.Ox-3 

1. Ox- 4 

1 . Ox-26 

31.7x-5 

30. Ox-5 

32.4x-5 

O.0x-5 

O.Qx-5 

0.73x-5 

14. 0/4 

11:40 

4,0 

8.0 

6.0 

0.0 

0.0 

10.0 

1.9x-6 

8 . Qx- 3 

1. Ox-4 

l.Ox-26 

31.7X-5 

30. Ox- 5 

32.4x- 5 

0 . 8x- 5 

O.Gx-5 

0.73X-5 

15. 0/4 

11:46 

4.0 

6.0 

6.0 

0.0 

0.0 

10.0 

2.5x-6 

B.Ox-3 

1.0x-4 

1. Ox-26 

31.7X-5 

30. Ox- 5 

32.4X-5 

O.Ox-5 

O.Ox-S 

0.73X-5 



